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We have developed a quantitative theory of Cooper 
pair pumping in gated one- dimensional arrays of 
Josephson junctions. The pumping accuracy is lim- 
ited by quantum tunneling of Cooper pairs out of 
the propagating potential well and by direct super- 
current flow through the array. Both corrections de- 
crease exponentially with the number N of junctions 
in the array, but give a serious limitation of accuracy 
for any practical array. The supercurrent at resonant 
gate voltages decreases with N only as sm(ip/N)/N, 
where (p is the Josephson phase difference across the 
array. 

PACS numbers: 73.23.Hk, 74.50.+r 

When a potential well propagates adiabatically 
along an electron system that is effectively one- 
dimensional, it carries with it additional electron den- 
sity and, in this way, induces a flow of dc electric 
current through the system. Such a pumping effect 
is observed in a large variety of mesoscopic systems 
ranging from small metallic tunnel junctions in the 
Coulomb blockade regime |Q, @, E|, to semiconductor 
quantum dots [||, and to one-dimensional ballistic 
channels §]. Th e propagation of potential well in 
these systems is arranged either directly through the 
propagation of real acoustoelectric wave [|[ || or by a 
few (two or more) phase-shifted gate voltages ||, || . 
Particular interest is attracted to the pumping regime 
when the potential well carries a quantized number 
m of electrons so that the induced current / is re- 
lated to the frequency /, with which the well crosses 
the system, by fundamental relation / = mef. The 
quantization of m is caused by the existence of the 
energy gap in the spectrum of the underlying electron 
system, which makes it possible for the potential well 



to be deep enough for precisely m electrons. Such an 
energy gap can be created either by the Coulomb in- 
teraction, as, for instance, in the Coulomb blockade 
pumps [0 ||, ||, or can be simply due to the discrete 
nature of single-particle states inside the well q, 0. 
In the case of Coulomb blockade pumps, precision of 
the pumped charge quantization is reaching the level 
sufficient for metrological applications [||. Different 
sources of inaccuracy in the pumps have been dis- 
cussed in §, | £| 0, | . 

Until presently, the pumping effect was studied 
mostly in normal systems, where transport is due to 
individual electrons. A timely motivation for study- 
ing Cooper pair transfer comes from quantum com- 
putation, where pumping can play an important role 
as an essential element of dynamics of quantum logic 
gates Jl2] . The aim of this work is to develop a quan- 
titative theory of pumping of Cooper pairs in one- 
dimensional arrays of superconducting tunnel junc- 
tions. In particular, we find fundamental corrections 
to the quantized pumping regime and show that they 
are unexpectedly large in arrays with a small num- 
ber of junctions. These large quantum corrections 
can also explain the fact that the first (and so far the 
only) experiment with pumping of Cooper pairs 
has failed to demonstrate accurate pumping. 

First we derive the general expression for the 
charge transferred through an array of N supercon- 
ducting tunnel junctions in the Coulomb blockade 
regime by adiabatic pumping of Cooper pairs. In 
the standard model such arrays are characterized by 
two energies, their charging energy He as a system 
of capacitors, and their energy associated with tun- 
neling . We assume that the characteristic energy 
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Eq of a Cooper pair in the array and the temperature 
(k^T) are both much lower than the superconducting 
energy gap of the electrodes. The first one replaces 
the usual condition in the normal Coulomb block- 
ade where the junction resistances should be larger 
than the quantum resistance. With these conditions 
fulfilled quasiparticlc tunneling is exponentially sup- 
pressed, while the tunneling energy of Cooper pairs 
in junction i reduces to a constant tunneling ampli- 
tude Eji/2. (Eji is also called Josephson coupling 
energy.) In this work, the bias voltage is set to be 
zero, and thus a constant Josephson phase difference 
ip is fixed across the array. We can then treat the 
two external electrodes of the array as one, so that 
effectively the array forms a loop and ip plays the role 
of external flux threading it. Then, the Hamiltonian 
of the N-pump is: 



N 



H = H c (n-q)^£ 

fc=l 



Eri: <\n){n + 5 k \e^l N +h.c 



(1) 



Here n = {ni, n-i, ■ ■ ■ , njv-i} arid q = 
{qi, qi, ■ ■ ■ , <Ziv-i} represent, respectively, the 
number of Cooper pairs on each island of the 
array, and qi, normalized by 2e, the charge injected 
to each island by the gate voltage V S i (Fig. 1). The 
term S k describes the change of n due to tunneling 
of one Cooper pair in the fcth junction. We will also 
need the operator of the current in the fcth junction: 



h = 



ieE 



2Ti 
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There are two mechanisms of Cooper pair trans- 
port in the array. One is the direct supercurrent 
through the whole array, another is pumping, the 
charge transfer in response to adiabatic variation of 
the injected charges q. To derive the general expres- 
sion for the total charge Q transferred during one 
pumping period, we introduce the basis of instan- 
taneous eigenstates |m) of the array for a given q. 
These are exact eigenstates if q is stationary. When 
q varies, there is a correction \8m) to state |m) of the 
form: 

\- \W\V q m) 



\5m) — ih 



(3) 



where ei >m are the eigenenergies of states \l, m). Thus 
the average current (Ik) in the fcth junction is differ- 
ent from the dc supercurrent (m\Ik\m) in the state 
\m), (Ik) = (m\Ik\m) + 2Re(m\Ik\5m) . Integrating 
over the pumping period t we obtain the charge Q, 
in units of 2e, transported through the array: 



Q= Ti~^J rfte '»(9W) + 2Re <j)(m\Q k \dr 



(4) 
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Figure 1: (a) A schematic drawing of a gated Joseph- 
son array of iV junctions. In pumping Cooper pairs gate 
voltages Vgi are operated cyclically, d are the capaci- 
tances of the junctions, and C g ; are gate capacitances. In 
a uniform pump C\ = C for all i = 1, 2, ... ,N. (b) A 
train of gate voltages to carry a charge in a pump. Here 
qi = CgiV^i/2e. 



Here we used the standard thermodynamic argument 
that the supercurrent in the state \m) can be ex- 
pressed as (2e/h)de m /d<p, and defined the operator 
of the normalized charge transferred in the fcth junc- 
tion, Q k = (l/2e) / dtl k (t): 



(m\Q k \l) 



ih (m\I k \l) 

2e El - E m 



(5) 



Diagonal matrix elements of Q k are not defined by 
eq. (||), but they do not contribute to Q. 

Few remarks should be made concerning eq. (|]). 
Since the state of the array at the end of the pumping 
cycle is exactly the same as in the beginning, charge 
conservation implies that Q is independent of the in- 
dex fc of the junction which is used in calculating Q. 
It is, however, a function of the array state \m). We 
did not include energy relaxation in the model and 
pumping is in principle possible even when \m) is an 
excited state. Below we consider only a more typi- 
cal situation when the array is in equilibrium so that 
at low temperatures \m) can only be a ground state. 
The last remark is that eq. ([|) shows that there is 
a close connection between the transferred charge Q 
and quantum mechanical phase of the state |m) accu- 
mulated during the cycle. Supercurrent contribution 
to Q is directly related to the dynamic part of the 
phase, Jq dt£ m (q(t))/Ti, while the pumped charge is 
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associated with the Berry's phase £ = i § (m\dm) |1E| . 

In the following quantitative analysis, the array is 
assumed to be uniform. The charging energy of the 
array is then: 



Ec 
' N 
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(6) 

where Ec = (2e) 2 /2C, C is the common capacitance 
of each junction in the array, and Uk = iik — qk- In 
the regime of accurate pumping the main contribu- 
tion to Q comes from the second term in (Q) while 
the supercurrent gives only small corrections limit- 
ing the pumping accuracy. The necessary condition 
for this regime to exist is Ej <C Ec, which we as- 
sume from now on. As in the case of Berry's phase, 
the second term in eq. ([|) does not vanish because the 
integration contour encloses the singularity where the 
energies of several charge states coincide. This degen- 
eracy occurs when — 1/N for all k. For such q, the 
array dynamics reduces to that of a particle on the 
N sites with equal energies forming a loop. The N 
eigenstates of such a particle are plane waves with en- 
ergies e fc = -Ejcos[((p-2irk)/N], k = 0,1, . ..,N-1, 
and the equilibrium supercurrent I (if) through the 
array is: 



Ic . f 

N Sm N 



(7) 



Here J c = 2eEj/h is the critical current of one junc- 
tion. Relation (^) should be continued periodically 
in ip beyond the interval [— 7r,7r], and I (if) exhibits 
cusps at if = ±7r. It also shows that the supercurrent 
decreases only as iV -2 at large N. 

Large supercurrent (Q) at resonant q means that 
the trajectory in g-space for accurate pumping should 
circle the degeneracy point sufficiently far away from 
it. It then successively brings in resonance the pairs 
of states that correspond to a Cooper pair occupy- 
ing two neighboring islands of the array (as illus- 
trated in Fig. 2 (a) for N = 3). If this process is 
slow, the Cooper pair is transported adiabatically be- 
tween the islands by the usual two-state level-crossing 
transitions that shift it along the array following the 
gate voltages. One Cooper pair is then transported 
through the array per cycle corresponding to a q- 
space trajectory circling once around the degeneracy 
point. One condition necessary for accurate pumping 
is that the probability of the Landau-Zener transi- 
tions to the excited states is negligible and the array 
remains in the minimum-energy state throughout the 
cycle. This condition limits the rate of pumping, 1/r, 
by the Josephson coupling energy, %/t -C Ej. How- 
ever, even then, i.e. in the regime of the present work, 



the pumping is not accurate due to the nonvanishing 
Ej/Ec- 

The gate voltages in an iV-pump are typically 
0, H sequences of triangular pulses shown in Fig. 
1 (b). The array is then completely translationally- 
invariant. The pumping cycle with triangular gates 
can be divided into N steps with Cooper pair trans- 
ported through one junction at each step. Instead of 
calculating the charge pumped in one junction dur- 
ing the whole cycle, Qp, we can equivalently sum up 
the charges Qj transferred in all junctions during one 
step of pumping: 

Qp = 2Re j>(m\Q k \dm) = EjLiQi ■ ( 8 ) 

Introducing variations 6(rij) of the average number 
of Cooper pairs on each island during one step, we 
can express the charge conservation (which follows 
from the Hamiltonian (Q) as S(rij) = Qj — Qj+i. 
Translational invariance implies that the distribution 
(rij) is just shifted by one island during the step, 
S(rij) = (rij-i) — (rij). For triangular gates we have 
precisely one Cooper pair in the array, EjL±(iij) = 1, 
and we then obtain: Ttf =1 Qj = 1 + N(Qjy — (tin))- 
Here the index N and the adjacent junction are arbi- 
trary. If we chose N such that the island is furthest 
away from the junction where the Cooper pair is tun- 
neling, we can neglect the occupation of the island in 
this equation and obtain: 



Qp = 1 + NQ N . 



(9) 



Now we can use eqs. (0) and (||) to calculate 
the charge Qn transferred during one pumping step 
through the junction that is most distant from the 
tunneling Cooper pair. Because of the energy differ- 
ence in the denominator of eq. (||), the main contri- 
bution to Qp arises at the resonances when a Cooper 
pair is transferred between the two islands. In this 
situation, there are two lowest-energy states with en- 
ergy separation on the order of Ej, and we can keep 
only the matrix elements of current between these 
two states in eq. (||). At the two islands, where the 
Cooper pair is transferred, the resonant states are 
given by the usual expressions of the two state sys- 
tems. To get a non-vanishing matrix element of the 
current in the A^th junction, we must "extend" the 
wavefunctions of the resonant states from the two is- 
lands occupied by the Cooper pair to this junction. 
We obtain by perturbation theory in Ej and by eqs. 
(|) and (|): 
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(10) 
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Thus the probability of Cooper pair tunneling lim- 
iting the pumping accuracy decreases roughly as 
(Ej I Eq) n ~ 2 with increasing N. It is physically clear 
that this conclusion should remain valid for non- 
uniform arrays also. Results of both numerical (from 
eqs. (|]) and (||)) and perturbative (eq. (|Io|)) calcu- 
lations are shown in Table I. 




Figure 2: (a) The states with minimum charging en- 
ergy of the uniform N = 3 pump on the (91,52) plane. 
The vector (711,712) denotes the stable configuration in- 
side each hexagon. A circular path with radius 5, and 
the triangular path of Fig. 1 (b) are shown. Only states 
contributing to the inaccuracy in charge transport in the 
leading order have been shown, (b) and (c) Numerical 
results of quantum inaccuracy of a uniform 3-pump for 
different values of Ej/Ec- The analytical result of eq. 
([n]), exact in the limit of small Ej/Ec, is shown by a 
solid line in (c). 

For N — 3, triangular gate voltages correspond 
to the triangular pumping trajectory on the (91,92) 
plane shown in Fig. 2. Another pumping method in 
the N = 3 pump jjl |l3| uses harmonic gate volt- 
ages, corresponding to a circular trajectory around 
the degeneracy point qi = q-x = 1/3. It is possible 
to calculate the pumped charge directly from eqs. (0) 



Table 1: Comparison between the numerically (top) and 
perturbatively (bottom) obtained transport inaccuracies 
(ip — 0) for uniform Cooper pair pumps with N = 3, 5, 
or 7 junctions and few values of Ej/Ec- 
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and (||). For tp = we obtain: 



Qr = l- 2 



-3V2S ±8 



1- 



where 5 = [( qi - 1/3) 2 + (q 2 - 1/3) 2 ] 1/2 is the ra- 
dius of the trajectory. The results of eqs. ( ]TT| ) and 
( JlC| ) for N = 3 almost coincide for the optimum ra- 
dius of 6 ~ 0.3. It should be noted that the quan- 
tum inaccuracy in pumping is very significant: it is 
more than 20 % at Ej/Eq — 0.03, which is a very 
small value. (Practically, Ej is limited from below by 
temperature, while maximum Eq is limited by min- 
imum feature size in fabrication.) The accurate co- 
herent pumping is thus practically impossible in the 
N = '& pumps. Figure 2 shows Qp calculated numer- 
ically from eq. (4) for ip = (no direct supercurrent 
present) as a function of 8. For small radii the charge 
is quadratic in 6, Qp/tt5 2 — (8i?c/27£ , j) 2 , as can be 
derived from eq. (0). At large 6's the pumped charge 
in Fig. 2 starts to decrease since the trajectory ap- 
proaches another degeneracy point at qi = q 2 = 2/3. 

Finally, we consider supercurrent in the regime of 
accurate pumping. In this regime, current is largest 
during the Cooper pair transition between a pair of 
islands of the array. The effective tunneling ampli- 
tude of this transition consists of two parts: the direct 
transition between the islands through one junction 
separating them, and the transition through the rest 
iV — 1 junctions. Interference between these two pro- 
cesses determines the phase-dependent part of the 
Cooper pair energy. We obtain the supercurrent of 
(N — Injunction tunneling by standard perturbation 
theory: 



Ej_ 
(11) 



1(9) = h 



E, 



N 



N-2 



(e 2 + Ejyn 



(JV-l) / Ej \ 
2(JV-2)! \2E C ) 



N-2 

sinyj, 



(12) 

where e is the energy difference between the Cooper 
pair states on the two islands of resonant transition. 
Current Jl2]) at resonance scales as (-Ej/ Eq) n ~ 2 , but 
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since the transition region represents only a fraction 
of the pumping cycle on the order of Ej / Eq , the over- 
all supercurrent contribution to the pumped charge 
(Q) scales as (Ej/Ec) N ^ 1 - This means that the inac- 
curacy is dominated by quantum tunneling of Cooper 
pairs (given by eqs. ( |lfj| ) or (|ll|)) even with nonzero 
ip if the pumping is not too slow. 

In conclusion, we have developed a quantitative 
theory of adiabatic Cooper pair transport in one- 
dimensional arrays of Josephson junctions. The the- 
ory predicts, among other things, that the quantum 
inaccuracy of the Cooper pair pumping in arrays with 
a small number of junctions is very large, the fact 
that can explain lack of success of the experimental 
attempt |l3| to pump Cooper pairs. 
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